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ABSTRACT 
A major hindrance to the use of algebraic techniques for the 
design of linear control systems is the non-linearity of the equations 
describing the system, The equations for a general class of control 
systems are shown to be linear in certain groupings of variables 
contained in, or added to, the systems, .A design technique based 
on the solution of these linear equations is outlined, and a 
digital computer program which will accomplish the design is pre- 
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Table of Symbols and Abbreviations 


s a complex variable (8 = o + Ae 

X(s) a transfer function; a ratio of polynomials in the 
variable 8s. 

Xn(s) the numerator of X(s) 

Xd(s) the denominator of X(s) 

Xn, the coefficient of the kth power of s in Xn(s) 

£(x,y, Z) a function of the variables x, y and z 

Kt variable forward gain 

T(s) forward transfer function 

K, variable feedback gain 

H(s) feedback transfer function 

C(s) filter compensation transfer function 

Ess Steady-state error 

G ob gain in decibels specified in a bandwidth requirement 

W frequency specified in a bandwidth requirement 

6 damping ratio, cosine of an angle in the s-plane 

oy undamped natural frequency, a radial distance in the 
s-plane 
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1. Introduction 

Linear control system design may be generally described as the 
synthesis of physical components or systems, active or passive, whose 
dynamic response will satisfy specified criteria. Present design 
procedure is essentially mathematical, since the response of dynamic 
systems is conveniently described by a set of differential equations. 
The problem is then to: 

1) formulate differential equations whose solution will 

describe the desired response, 

2) relate coefficients of the equations to the parameters of 

a physical system, 

3) obtain a solution of the differential equations in terms 

of these parameters. 
Production of systems designed in this manner is frequently economi- 
cally unfeasible, since components required may be expensive to man- 
ufacture in limited numbers. A more practical approach is to utilize 
readily available components, combining them and adjusting their 
variable parameters to obtain a system with the desired dynamic 
response, 

Analysis of systems thus synthesized may be accomplished in the 
time domain by direct solution of the system differential equations, 
but the more usual approach is to perform analysis in the complex 
frequency domain (s-plane), by means of the Laplace Transformation. 
System equations thus transformed may be approached with several well- 
known analysis techniques, notably the Evan's root locus procedure, 
or in the frequency domain (s = 4%)» the Bode, Nyquist or Nichols 


methods. Application of these techniques to the design procedure is 
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restricted, however, since they are not easily described analytically, 
they are essentially graphical and therefore time-consuming, and the 
approach to design is largely trial and error, All suffer the further 
disadvantage that changes in system response may be predicted on the 
basis of variation of only one parameter. 

Several recent techniques, to be discussed in further detail, 
have succeeded in overcoming the more obvious difficulties listed 
above. The transformed differential equations describing the system 
and the specifications prescribing system response are reduced to a 
set of algebraic functions of the variable parameters, The algebraic 
relationships produced, however, are non-linear, and for complex 
Systems may prove difficult or impossible to solve. This paper out- 
lines a numerical procedure for the solution of these non-linear 
equations, and presents a digital computer program which obtains 


solutions for a general class of control device problems. 
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2. Algebraic Methods 

Application of algebraic techniques requires the reduction of 
the transformed system equations to a "transfer function," usually 
the ratio of two polynomials in the variable s, which relates the 
output of the system to an arbitrary input. The roots of the denom- 
inator polynomial, known as the "characteristic equation," determine 
the stability and response time of the system, factors of principal 
concern to the designer. For the system to be stable, i.e. finite 
response to a finite input, it is necessary that all roots of the 
characteristic equation have negative real parts. It can be shown 
that the time response of the system is determined by the values of 
the roots, and is frequently dominated by the root or pair of roots 
whose real parts have the least negative value, i.e. when plotted in 
the s-plane, the roots located closest to the origin. The dominant 
roots are designated in polar coordinates by their radial distance 
from the origin (wy), and the cosine (¢) of the angle formed by the 


second quadrant radius vector and the negative real axis. 


¢ UW 


s-plane 


Ss @ oF .w 
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Figure 1 
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Method of Mitrovic' 


The characteristic equation is a polynomial of the form 


n nel 2 1 
F(s) Fas + F i-l 8 + 2--= + F,8 +- Fi8 a F (2.1) 
or in summation notation 
a kk 
F(s) = 2 Eis® (229 
k=0 


where the polynomial coefficients (F,.'8) are functions of the system 
parameters, Mitrovic'[1] has shown, through a transform of variable 
and s-plane mapping procedure, that two coefficients of F(s) may be 
Stated as functions of the remaining coefficients and of the — and 
Wn of a pair of specified roots. The result of this procedure is two 


equations of the form 


k 
E (+1)" Fun O47) = 0 (2.3) 


y (-1) Fy Wn (8) = 0 (2.4) 
k=0 
where oy (E) and o,, 8) are recurring Chebishev functions, which for 


a specified — may be precomputed from the relationship 
b14 66) 726, CE) + 4,.,) = 0 (25) 


with ¢ ©) = 0 
and >, ) = «] 
If numerical values are assigned to all but two of the coeffi- 


n? equations (2.3) and (2.4) are reduced 


cients (FL) and to € and w 
to a pair of simultaneous equations in two unknown variables, i.e. 
the coefficients which have not been specified. Mitrovic' originally 


chose the coefficients E and F. to be variables, arriving at the 


1 


relations 
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F = f, (Fy Ey ui) (2.6) 


Oo 
Fy ms f, (FL Ss wi) (2.7) 


Later, Sikjak generalized the Mitrovic' method to show that a solution 
can be obtained for any pair of coefficients [2]. Thus, a relationship 
between the coefficients of the characteristic equation, which are 
functions of the system parameters, and the location of a desired pair 
of roots is established, Designating the variable system parameters 


as (a, > Ay, ooo a)» we have for any two coefficients 


F 
x 


£,(a,> ao» eet a.) © £, (FL > E wn) (2.8) 


Fy — £5 (a,> ay» a a.) = £6 (FL » Ey» Wn) (232) 


The process may be repeated to specify additional pairs of roots, 
thus establishing additional relationships of the form of (2.8) and 
(2.9). If, however, not all roots of the characteristic equation 
are specified, the method provides no information on the location of 
undesignated roots, and does not guarantee either stability or the 
dominance of a desired pair of roots. 
Coefficient Matching 
If all of the roots of F(s) are designated, the characteristic 


equation may be stated as 


n 
F(s) = « (8+ r) (2.10) 
k=1 
where (r)> Tos 797 r) are the roots of F(s). The resultant coeffi- 


cients of the expanded polynomial are functions of the roots. 


n-1, ae 


n 
F(s) = £ (8) >8os7r or) 8 + fel (r)»Too77 "x )8 


Z 
-+ £, (44985 o7778)8 +£, (4, Tose 9t D8 + f(r, »toreets) 
(2.11) 
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where £, (ry s¥o> --- x) is the sum of the products of all possible 


combinations of (i> pote} taken k at a time. 


9° 
For example: 
FY (s) = f. (X2T5> ==) 
ie at 7 eT aan cS A a la 9, 
+ (ry 84r,---T 1) + ---- (2.12) 


In summation notation 


x n 
Fe) wt [™k (r5)] _x Sent (2 3) 
i=l jai k! (mek)? 


Since all roots are designated, the F,'s are fixed constants. The 


k 


result is two expressions for F(s) 


n 


F(s) = ¢f £,(a,,a,-¢ra )st Parameters (2.14) 
p Rite, Laney Z 
i=0 
. i 
F(s) = ¢: g.(%,,r,---r_)8 Roots (2.15) 
i=0 2 Loz n 


In order that these expressions be identically equal, coefficients 


of the same powers of 8 must be equal. 


- .. k = k 
F,8 = £, (a, 2427-748 = (ry sFos7--8 8 (2.16) 


The resultant system is composed of n simultaneous equations in z un- 
knowns. Complete solution of the system of equations requires that 
the number of unknown quantities (z), the variable parameters, be equal 
to the number of equations (n), the order of the characteristic equa- 
tion. If insufficient variables are available, additional parameters 
must be inserted in the problem in order that this requirement be 
satisfied. 

The principal advantage of designation of all roots is the guar- 


antee of stability and dominance of desired roots. By designating a 


16 


pair of roots whose values give the desired time response, and loca- 
ting the remaining roots in such a way as to inSure dominance of the 
original pair and stability of the system, the designer completely 
specifies the time response of the system. The principal disadvant- 
age is the necessity for addition of arbitrary parameters sometimes 
required to permit solution of the simultaneous equations. 

A disadvantage common to all algebraic methods is the inherent 
non-linearity of the resultant equations. The expressions usually 
contain products and/or powers of the variable parameters, and 
assistance of a digital computer is required for solution of all 


but the most simple cases. 
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3. Algebraic Solution of a General Problem 

Given a dynamic device whose transformed time domain equation 
may be expressed as the product of a variable gain (Kt) and a trans- 
fer function T(s), and a feedback device expressible as a variable 
feedback gain (Kf) and the transfer function H(s), the transfer 
function F(s) of the "closed loop" control system comprising the two 


components, as arranged in Figure 2, is 


Output (s) .. 2 Kt T(s) 
Tnpeeece) F(s) 1 + Kt T(s) Kf H(s) (3.1) 






Output (s) 





Figure 2 


Basic Feedback Control System 


Both T(s) and H(s) are the ratio of polynomials in S. 





a aol 2 
8s > @ 2 @ @ 
Tn(s) In, ft fn 41* + +: Tn,8 + hil 15 Tay 
oe Td(s) - Td 3” + Td gb-1 + eooe + Td 3? + Td.s + Td 
b b-l 2 ? re) 
s° (wn.e* + nis + Ha) (3.2) 
Hn(s Mo vol = : 
(OS ant. Won BD 
Hd.s + Hd.s + Hd 
2 1] re) 
Kt Tn(s) 
and F(s) = Td(s) (3.3) 


1 + Kt Tn(s) Kf Hn(s) 
Td(s) Hd(s) 


F(s) reduces to the ratio of two polynomials 
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> Sree Kt Tn(s)°Hd(s 
Se) race) Td(s)-Hd(s) + Kt K£ Tn(s)-Hn(s) (3.4) 


Forming the product of polynomials where indicated defines: 


Tn Hd(s) Tn(s)*Hd(s) 


Td Hd(s) = #£Td(s)-Hd(s) 


Tn Hn(s) Tn(s)°Hn(s) 
The characteristic equation of the system - Fd(s) is the sum of the 
polynomials TdHd(s) and KtKfTnHn(s). 


Fd(s) 


TdHd(s) + KtKfTnHn(s) 


TdHd s” a5 TdHd_ Et + ---- + TdHd.s + TdHd | (3.5) 


She 1 
y-l 


y 
+ KtK£(TnHn + TnH 
( y® ayes 


+ «ae + TnHn,s + TaHn , 


In most devices used for control systems, particularly higher order 
systems, the order of the numerator polynomial is much less than the 
order of the denominator. Feedback devices seldom have numerator 
order more than two greater than denominator order, i.e. second 
derivative feedback. Thus the order of the characteristic equation, 
and the number of roots to be designated, is determined by the order 
of the TdHd(s) product. By expanding the designated roots into a 
polynomial R(s), and equating coefficients, the system of simultaneous 


equations for Fd(s) becomes 


x . x 
(TdHd_ ) 8 = R, é 
x=] x=] 
= 8 
(TdHd_,) 8 Raa 
(TdHd, = + KtKET niin, yay a RS” 
sy~1 = gy-l 
(TdHdy_, + KeKéTniny_ 1) om 
: 2 ” 2 
(TdHd, + KtKfTnin, )8 = R,8 
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(TdHd, = KtKffnkn, )s = R.8 


(TdHd 
Oo 


i) 


KtKfinlin, ) = R (3.6) 


If all polynomials are rationalized, 

TdHd | = R. = 1 
and the system consists of x equations in the two unknown parameters 
Kf and Kt. In order to obtain an algebraic solution, an additional 
(x-2) variables must be inserted in the system. The addition of 
"cascade compensation" in the form of a multi-section filter net- 
work permits introduction of the necessary mathematical variables 
and is also an economical method of inserting variable parameters 


in the physical plant. 


Cascade Compensation 


Input(s) 






Output (s) 


Figure 3 


Basic Feedback Control System with Cascade Compensation 
The transformed output/input equations C(s) of a resistance- 
capacitance filter network are expressible as the ratio of two 
equal order polynomials, the order being equal to the number of 


filter sections employed, 


af 


m 
Cds oa Cd 4° 


Cn —— Cn + wae + Cana” +: Cn,s + Cn 
m m a 1 








C(s) = ae —— j (3.7) 
* ab eae + Cd,s° + Cd,8 + Cd. 
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The transfer function of the compensated closed loop system, for 


the configuration shown in Figure 3, becomes 


Output (s) _ ns Kt T(s)eC(s) 
Input (s) RK) 1 + Kt T(s)eC(s) KE H(s) (3.8) 


which reduces to 


* Fn(s) - Kt TnHd(s)e°Cn(s) 
ate) Fd(s) TdHd(s)*Cd(s) + KtKf£ TnHn(s)eCn(s) (3.9) 


The order of the characteristic equation is now determined by the 
TdHd(s)eCd(s) product, thus the addition of each filter stage 
increases by one; the order of Fd(s), the number of simultaneous 
equations to be solved, and the number of independent variables 
required. If the order of TdHd(s) is (x), and (m) filter sections 
are used, the order of Fd(s) will be (x +m). The addition of each 
filter section (i) introduces two independent variables into the 
systen, Cn, and Cd,- Thus (m) filter sections introduce (2m) 
variables. 

For solution of the system of equations, the product of the 
independent parameters (KtKf£) may be considered as a single 
algebraic variable. Addition of (m) filter sections to a system 
whose TdHd(s) is of order (x), produces (x + m) equations in 
(2m + 1) variables, Equating the number of variables and number 
of equations 

x + m = 2m +1 
me x-l (3.10) 
For the system given in Figure 3, the number of filter sections 
required is determined by the product of the denominator polynomials 


of the forward and feedback path transfer functions. 
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The system of simultaneous equations for the compensated system 
is obtained by designating the required (x + m) roots, expanding the 
roots to form the polynomial R(s), forming the TdHd(s)*Cd(s) and 
KtKfTnHn(s)*Cn(s) products, and equating coefficients. Table 1 
illustrates, in matrix form, the equations for a system having a 
third order Td(s)+Hd(s) product, a first order Tn(s)*Hn(s) product, 
and the required second order filter. The format of the matrix 
equations is the same for higher order systems if the requirements 
that 

1) The number of filter sections used be one less than the 

order of Td(s)-Hd(s), 

and 2) The order of Tn(s)*Hn(s) be less than the order of 
Td(s)+Hd(s) 

are fulfilled. 

Note that the equations for systems of this type are linear in 
certain individual parameters (Cd; Cd, one Cd); and in groupings 


of the other parameters (KtkfCn, KtKfCn,, --- KtkfCn_). Thus the 


A 
system of equations may be solved by linear methods for these 
individual and grouped parameters, 


Again, if all polynomials are considered rationalized, 


TdHd Cd. = R in = 1 


so Cd =1. If Cn is defined as equal to 1, the group variable 
KtKfCn becomes simply the product KtKf. Having solved the linear 
system and obtained the coefficients of the denominator polynomial 


(Cd; Cd,, --- Cd 1) of C(s) directly, the coefficients of the 


]? 
numerator polynomial may be obtained by dividing the grouped variables 


(KtKfCn_, KtKfCn,, ---- KtKfCn_.,), by the value of KtKfCn , i.e. 
oO m=] m 


1? 
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by KtKf. With the numerator and denominator polynomials 
available, the filter design becomes a problem in passive network 
synthesis. 

A solution has thus been obtained for all physical parameters 
except the individual gains, Kt and Kf. An additional specification 
is required in order to obtain an independent solution for values 
ef the gains, The additional criteria considered are steady state 
error and bandwidth. 

Steady - State Error (Ess) 

The error of a system is defined as the difference between 
an input command signal and the system response, 

Error (s) = Input (s) - Output (s) (3.11) 
From Equation 3.8 


Output (s) = Kt T(s)*C(s 
I+ Kt T(s)°C(s) KE H(s) repute (s) 


The error as a function of the input signal becomes 


1 + Kt T(s)*C(s) {KE H(s 


a 
1 + Kt T(s)-Cfs) Kf H(s) Input (8) (3.12) 


Error (s) 


which reduces to 


Error (s) _ _1 Cd(s) + (s) = KtTnHd(s)-Cn(s) 


Input (s) Tdd(s)eCd@s) + KtKf Tnlin(s)=cn(s) (ate) 





For a unit step displacement input, Input (s) is l/s. Error (s) 


for a step input becomes 


1 Kt TnHd(s)*Cn(s 
Ermas (3) hag E ” ‘TdHd(s)eCd(s) + KtKE aa OES) (3.14) 


Steady-state error is defined as the error of the system after all 
transient effects have taken place, or mathematically, as time 
expired after application of the input signal approaches infinity. 


Application of the Final Value Theorem to Equation 3.14 gives 
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Ess = lim s Error (s) 


s>o 
Ess = lim - Kt TnHd(s)*Cn(s) 
SO TdHd(s)eCd(s) + KtK£ TnHn(s) «Cn(s) 


TdHd «Cd + TnHn KtKfCn - TnHd »« KtCn 
re) re) re) O re) oO 


TdHd «Cd + TnHne ktkfCn (3.15) 
Oo Oo Oo O 


Note that steady-state error as expressed above is defined for unity 
input. Error is relative to input, and maximum steady-state error 
desired in the system is specified as a fraction of the input signal. 

Steady-state error is found to be a function of two variables 
previously defined, Cd. and KtKfCn |, and an additional variable, 
KtCn.- Values for Cds KtK£Cn | and KtKf having been previously 
obtained, the solution for Kt Cn, is 


(1 - Ess) [TdHd »Cd_ + TnHn~-KtKfCn J 
Oo Oo Oo Oo 





Kt Cn, = Tad (3.16) 
Oo 
KtKf£Cn 
te) KtKf£ 
and Kf “ttn, ; Kt KE (3.17) 


Values of all unknown parameters have been obtained, and a system 
is synthesized which will meet specified requirements for stability, 
root dominance, and steady-state error. 

Bandwidth (BW) 

The bandwidth of a system is defined as the frequency range 
over which the response of a control system is within acceptable 
limits [1]. It is frequently given as that frequency (w) at which 
the magnitudes of the output response and the input signal have a 
specified ratio (G), usually expressed in decibels (db). 

If F(s) is the closed loop transfer function of a control 


system, bandwidth may be found from 
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output (. w) 
Cap 7 20 10819 \ Guput Go) emeiee: Weil (,~)] (3.18) 


where 4 is the desired bandwidth, 


Fn(,u) Fn(,u)* + === + Fn, (,u) + Fn, (,4) + Fn, 


cl 5, Fa.) 


b 2 
Fa, (ju) + me + Fd, (,0)"+ Fay (,4) + Fe, 


16; 


a 2 Fn(y)real_ + j Fn(,)imag, 
Fd(w)real + j Fd(w)imag. 


Fn(w) real ‘ + (Fn(w)imag. 27) 1/2 

G = 9 9 (3989) 
(Fd(w)real)” + (Fd(~) imag.) 

The bandwidth expression is seen to be highly non-linear, however, 

this obstacle is not insurmountable. Recalling that 


Kt TnHd(,)-Cn(, ¥) 


nw) os 
TdHd ( , w)- Cd, w) + TaHin( , w)e ReK£Cn(, 4) 


and that values for the variables Cn,» Cd. and KeKiCn are available 
from the solution of the simultaneous equations, Kt may be obtained 
from the bandwidth criterion. If gain (G) is specified at frequency 
(4), 


TdHd(.w)eCd(, 4) + TnHn(, “)eKtK£Cn(. %) 
Kt = G Eg ee et 


Tad ( <u) eCn¢, u) (3.20) 


and K£ = KtKet (3.21) 
Kt 
Again, values of all unknown parameters are obtained, and the system 
design will meet stability, root dominance and bandwidth specifications. 
Multiple Specifications 
Either a steady-state error or a bandwidth specification is 
sufficient in itself to permit solution for all variable parameters 


and design of the system. However, if both error and bandwidth 
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are specified, solution for one does not guarantee satisfaction of 
the other requirement, A slightly different approach must be taken. 
No requirement has been placed on the location of the non- 
dominant, or non-specified roots, other than that they have negative 
real parts, and‘ that they be sufficiently removed from the origin of 
the s-plane to insure dominance of the desired roots. This gives the 
designer a great deal of latitude in the choice of extraneous roots 
which are added to the problem as a result of cascade compensation. 
The system of linear equations used in the solution for para- 
meters is obtained by equating coefficients of the characteristic 
equation. On one side of each equation is an algebraic function 
of the variable parameters; on the other side is the numerical 
coefficient found by expanding the roots to form the polynomial 
R(s). Changing the location of one or more roots changes the values 
of the variable parameters, Parameter values obtained are those 
which will place the roots of the characteristic equation exactly 
where specified. The steady-state error specification is satisfied 
by determining the relationship of the forward gain (Kt) to other 
system parameters, The bandwidth is found from the equation 


Kt TnHd ( w) Cn. w) 
RC, 0) 


G (3,22) 


Thus the bandwidth of the system is a function of the location of 
the roots, but the steady-state error solution is not. By judicious 
choice of his root locations, the designer may be able to meet both 
error and bandwidth specifications, or to determine whether a band- 
width specification conflicts with a stability or root dominance 


criterion, 
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The process of finding the proper root locations is iterative, 
laborious, and time-consuming, requiring a complete solution of the 
system equations after each movement of the roots. For higher order 


systems, the assistance of a digital computer may be necessary. 
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4, Limitations of the Method 
Steady-State Error for a Ramp Input 
Steady-state error, as defined in Section 3, is the limit as 
S approaches zero of the product of s and the error function 


Error (s), where 


Error(s) = Input(s) TdHd(s)eCd(s) + KtKfTnHn(s)*Cn(s) - KtTnHd(s)+Cn | 


TdHd(s)*Cd(s) + KtKf TnHn(s)Cn(s) 
For a unit ramp input, i.e. constant velocity signal of unity 


magnitude, Input (s) = L/s*. Steady-state error becomes 


Fie alti TdHd(s)*Cd(s) -}- KtKfTnHn(s)¢Cn(s) * Kt Tala (s)*Cn(s) | 
S30 82 TdHd(s)eCd(s) + KtKf TnHn(s)Cn(s) (4.1) 


This limit approaches infinity as s approaches zero, unless an s can 
be factored from the numerator polynomial of the term in brackets, i.e. 


. TdHdeCd + TnHneKtKfCn - TnHd«KtCn = O (4.2) 
o 6 O Oo Oo Oo Oo 


Equation 4.2 may be set equal to zero in several ways: 

1) The possibility of requiring that Cn, or Cd. be zero must 
be disregarded, as this would fix the value of a variable 
required for solution of the system of simultaneous equa- 
tions. 

2) Hn, may be required to be zero, i.e. tachometer feedback, 
but the steady-state error expression resulting from this 
procedure introduces two new variables KtCn and KtCn, 
in a single equation 


Ess = £(Cd., Cd KtKxCn_., KtCn_, KtCn 


L » 


and solution for independent Kt is impossible. 


Kt Cn, Kt Cn, 
Kt men ae # Kt nel. a 
re) 1 
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3) Td | may be required to be zero, i.e. a pole at the origin 
of the open loop transfer function, Then 


Tn Hn KtKfCn -~ Tn Hd KtCn = Q 
o Oo oO o Oo oO 


and Hd s i.e. unity feedback. 


Kf = Oo 
Hn 
o 
If this procedure is used Kf is no longer an independent 
variable and solution is impossible, 
4) If Hd | is required to be zero, 


Tn Hn kKtKfCn = QO 
oO re) 


or Tn, must also be zero, and an s factors from the denomin- 
ator as well as from the numerator of Error (s). 
To date, the author has found no expression for steady-state 
error to a ramp input which is compatible with the system of 
simultaneous equations previously derived. 


Multi-path Feedback 


Input(s) 
Kt T(s) C(s) 
+ Output(s) 





Figure 4 
Control System with Two Feedback Paths 
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The transfer function for the control system illustrated in 


Figure 4 may be written 


Kt T(s)*C(s 
F(s) i+ Kt T(s)*C(s) (Kf H(s) + Ke} (4.3) 
which reduces to 
¥ Kt Tn Hd(s)eCn(s 
KS) S TdHd(s)*Cd(s) + KtK£ TnHn(s)eCn(s) + KtKg TnHd(s)eCn(s) (4.4) 


The system of simultaneous equations obtained by equating coefficients 
of the characteristic equation will be linear in the variables (Cd; 
Cd 


w-- Cd_), (KtK£Cn,, KtKk£Cn,, --- KtK£) and (KeKgCn,, KtKgCn,, --- 


1? Ve 
KtKg). The equations may be solved to obtain values for the 
variables composed of grouped parameters, However, for this systen, 
and other multi-path feedback systems of various configurations, it 
has not proved possible to sclve for independent values of the 
independent parameters... For example, in the system above KtKfCn 
KtKf, KtKgCn_, and KtKg are independent variables. But the value of 


Cn, obtained from 


Cn = KtKfCn 
re) ’e) 
KtKf£ 


is not necessarily the value of Cn) obtained from 


Cn, = KtKgCn 
KtKg 


Systems of this type have defied solution by this method, 
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5. Sample Solutions 


Example 5.1 


I(s) 1 
zt * (stl) (8-2) (848) (s+15) 


Kf 


Figure 5.1 
Block Diagram for Example 5.1 


The open loop transfer function given above is unstable due 
to the pole in the right half plane (s = 2). Design a filter 
compensator and find values for —— and feedback gain which 
will make the system shown in Figure 5.1 absolutely stable. 


Dominant roots are specified by € = 0.707 and o = 2a « 


No additional steady-state error or bandwidth requirements 


are specified, 
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~ 0(8) 


INPUT DATA 


RHEKRHH HAHAH 


CPER LOCP TRARBRBER MEU eign = 1 (5S) 


NUMERATOR CF OPEN LOOP 


TRANSFER FUNCTION - TNS) 


TN(S)= 1.°°RGE494 


BENGE itt OF CPREN -E@Gr 
TRANSFER FUNCTION - TODS) 


NUMBER GF ROCTS = 4 


ROOT NC. J Es 3 4 


Re APT eee +) S2UEt owe tol SL OE TO Z 


IMAG PT. e tia oe re renee oS eet > (: 


COEFFICTENTS OF FEXPANDEDMPCL YG =ieAt 


TOGS)= 1.57686 +20.Saee 


t 2.28800 + SS 


pms 6 SUC EO eS ( 


4-1] ,66060¢4+°02 St 1) 


a Aer EO a 


NOTE = PROGRAM HAS FP CAD PURES 


IN THE RIGHT HALF OF THe 
COMPLEX PLANE. 
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FEECBACK PATH TRANSFER FUNCTION - H(S) 
(0) OROER FEEDBACK IS SPECIFIED. 
NUMERATOR OF FEEDBACK PATH 
TRANSFER FUNCTION - HN(S) 
HN(S)= 1.0000E6+00 
DENOMINATOR OF FEEDBACK PATH 
TRANSFER FUNCTION - HD(S) 


HD(S)= 1.9000E+00 


CLOSEC LOOP TRANSFER FUNCTION 
G(S) - UNCOMPENSATED 


NUMERATOR CF CLOSED LOOP 
TRANSFER FUNCTION - GNCS) 


GN(S)= KT#® 1.0G00E+00 


CENOMINATOR CF CLOSED LOOP 
TRANSFER FUNCTION - GO(S) 


GD(S)= 1.90CGFFHCO SC 4) 
+ 2.20CGG701 Site 
+ 9.50C0L+01 St 2) 
+-~1.6600E+02 S( 1) 


+ (-2.4000E+02 + KT#KX# 1.0000E+4+00) 


ST EAD i ee ee Ser eC ERC AT LG 


NO STEADY STAT= ERROER REQUIREMENT 
IS SPECTELED. 


BANDWIDTH SPECIEICARIO 


NO BANDWIDTH RFOUIREMENT I[S SPECIFIED. 


NOTE —- DATA INOICATES “I0 BANDWIDTH OR PARAMETER 


DEPENDENT SIPADY SPATE ERROR SPEC TEICATIONS. 
UNDER THESE CONDITIONS TRE PROGRAM WILL NOT 
SGLVE FOR FORWARD GAIN AND FFEDBACK GAIN 
SEPARATELY. SU CUTION TS (CR THES PvE Torr 
THE GAINS. (KT#KX) 


RIEU See SPECIE] Ee amg 


CEUSEE EOP = sieuT [ON 


LEWN = «eer Clie OMEGA = 2.5090 Etue 


NUMBER OF ROOTS = 2 


ROOT NO. ] 2 


REALSPh.s =ellSetrl: -~.i SE+8l 


IMAG PT. SOGEtOL =. l Get ul 
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SOLUTION OF THE PROBLEM 


\PBRRAARSZRALAARAAALARLALLR ESL SS SE 


PROGRAM REQUIRES 7 ROOTS FOR SOLUTION. 


SOLUTION OF THE SYSTEM OF LINEAR EQUATIONS 


VAR TABLE PARAMETERS RELATIVE ERROR 

VARIABLE ( 1) CBO) 7.358E-09 PERCENT 
VARIABLE ( 2) CD(1) .003E-01 PERCENT 
VARIABLE ( 3) Co(2) ~O00E-01 PERCENT 
VARIABLE ( 4) CD(3) -4.180E-09 PERCENT 
VARIABLE ( 5) KT#KX#CN (0) -1.032E-07 PERCENT 
VARIABLE ( 6) KT#KX#CN(1) 3.466E-07 PERCENT 
VARIABLE ( 7) KT#KX#CN (2) 1.393E-06 PERCENT 
VARIABLE ( 8) KT#KX#CN (3) 2 a@le=06 PERCENT 


NOTE - RELATIVE ERROR STATED ABOVE 
MEASURES THE ACCURACY OF THE 
SOLUTION FOR EACH VARIABLE. 
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NOTE 


SOLUTION FOR INDIVIDUAL PARAMETERS 


NUMERATOR OF CASCADE COMPENSATION 


CNOS)= LeGacorer® $( 3) 


+ 1.1867£64+91 S( 2) 


+ 3.43520ETOT Sil 


+2. 7105Et01 


FERUS OF ait bk be Ges 


ROOT NO. l a 3 


REAC Pt. =SeOr+Cl =. Teer) Seeeesel 


IMAG PT. ac tle oN IF +2 wee 


- ACCURACY OF ROUT FINDING SUBRCUTINE 


HAS BEEN CHECKED BY EXPANDING THESE 
ROOTS, AND COMPARING THE COEFFICIENTS 
OBTAINED «id Tee THE COBEFAGIENTSeSE THE 
ORIGINAL POLYNOMIAL. 

ERROR NECTED IS LESS THAN CNE PERCENT 
FOR ALL COEFFICIENTS. SOLUTION FOR 
RCUTS IS CONSIDERED ACCURATE. 
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DENOMINATOR OF CASCADE COMPENSATION 
CO(S)= 1.0000E+00 S( 3) 

: 1.5117E+01 St 2) 

+ 1.4418E+02 St 1} 


+ 2.9161E+02 


POLES OF THE FILTER 
ROOT NO. 1 a 3 
REAL PT. -.63E+01 -.63E+01 -.26E+01 


IMAG PT. -.85E+01 .85€&401 .00E+00 
es ate he 
NOTE =! AGCURACY HAS BEEN CHECKEO,AS BEFORE. 
“ERROR NOTED IS LESS THAN ONE PERCENT 
FOR ALL COEFFICIENTS. SOLUTION FOR 
syaRQOTS IS CONSTDERED ACCURATE. 


PRODUCT OF GAINS (KT#KX) = 5.4508E+03 
7. oe 


FORWARD’ GAIN (KT) IS 5-4508E+03, 
WITH UNITY FEEDBACK GAIN. (KX=1.0) 
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CLOSED LOCP TRANSFER FUNCTION 
OF SOLVED PROBLEM 
NUMERATOR POLYNOMIAL 
GN(S)= 5.4568F+03 S¢ 3) 
+ 6.4686E+C4 S( 2) 
+) lob TOC oy Seay 


tale 4 TASES 


ZEROS OF THE TRANSFER FUNCTICN 

ROOT NO. 1 2 3 

REAL PT. =-.80+01 =SS ea —.25 tee 

IMAG PTI. silefel +0 ~UCE+O2 COL +59 
NOTE = ACCURACY FAS. BEEN CHECKED SAS SREEO RE. 

| ERROR NOTED TS LESS THAN TENE @eearcenT 


POUR ALL COEFFICIENTS Se SGL U TIGinerOR 
RCOTS IS CONSIDERED ACCURATE. 
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CENOMINATOR POLYNOMIAL 
GD(S)= LeCocor4dO St 7) 
+ 3.71186+61 S( 6) 
+ Sov iit 2S (5 | 
+ 4,.7338E+C3 S( 4) 
+ 2.281464+uU4 S( 3) 
+ 6.4826F+°4 S( 2) 
+ 1.94Q06F+05 S( 1) 
+ (.7761E+04% 
POLES OF TRE SYSTEM AND ROOTS 
CF THE CRKARACTERISTIC EQUATICN 
ROOT NO. l 2 3 4 5 
REAGwPis =sroctol -.PoerToL “—~sosetul -.53EFO]l —. (ieee 


IMAG PT. wR L -.-_)8Eae 1 elBEt99 -. LBEF0O e L8E+00 


ROOT NC. 6 T 
REAL PT. —-e7lEtOlL -. 886481 


IMAG PT . -.-18F+4209 SIOETI 
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Example 5,2 












I(s ee 
yee oa d 


8 + 338° + 3368- + 1880s + 7200 O(s) 





KE (gee 
ss’ + 30s + 216 


Figure 5.2 
Block Diagram for Example 5.2 


Cascade compensation and the tachometer feedback device given 
are to be used in the closed loop system shown in Figure 5.2, The 


maximum permissible steady-state error for a step input is to be 


10 percent. 


Dominant roots are specified by ~ = 0.5 and y= Si. 


No bandwidth requirement is specified, 
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INPUT DATA 


SH#SSeeeeaeaeae 


OPEN LOOP TRANSFER FUNCTION - T(S) 
NUMERATOR OF OPEN LOOP 
TRANSFER FUNCTION - TN(S) 
TIN(S)= 1.0000E+00 S( 1) 


+ 5.0000E+00 


DENOMINATOR OF OPEN LOOP 
TRANSFER FUNCTION - TOCS) 


TOCS)= 1.9000E€+00 St 4) 
+ 3.3000F+C1 S( 3) 
+ 3.3600F+92 S( 2) 
+ 1.8800F+03 S( 1) 


+ 7.2000t +03 


ROOFS OF TOS) 
ROOT NO. 1 2 3 & 
REAL PT.~ --90EtOL --20E+02 -.20E+01 -.20E+01 


IMAG PT. eOOE+00 eOOE+00 -.60E+01 e60E+01 
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FEECRACK PATH TRANSFER FUNCTION - H(S) 

(1) ORDER FEEDBACK IS. Ses olEmcon 
NUMERATOR UF FEEDBACK PATH 
TRANSFER FUNCTION - HN(S) 
HN(S)= 1.000045 S( 2) 

+ 7.CG90GOlP+CC "St 1} 

+ seg =") 
CENOMINATOR OF FEEOKACK PATH 
TRANSFER FUNCTIGN - HD(S) 
HD(S)= 1.40COR+70 St 2) 

7 cepeece+) | Se 


+ 2.160C0r+52 
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CLOSER LOOP TRANSFER FUNCTION 
GUS) =JUNCOMBENSATLED 


NUMERATOR CF CLOSEN LOOP 

TRANSFER FUNCTION - GN(S) 

GN(S)= KT#( 1.0007£+00 S( 3) 
+ 3.96 OE tals 2 ) 
+ 3.660CE+a2 S(t 1) 


+ Lies 60S es) 


CENOMINATOR CF CLOSED LOOP 
TRANSFER FUNCTION - GDCS) 


GDIS)= Lo yf ee oe) 
+ 6. 30600 G+ )eeSt 5) 
+ 1.5420F+%93 SC 4) 
+ ( 1.9U088E+14 + KT#KX@# 1. GQUDEt 7) SC 3) 
+ ( 1.3618F+°5 + KT#KX# 1.2000F+°1) SC 2) 
tunel He 22 PEt. 5 eee se Steet i ees 1) 


Hot 1 .5552E+ 6 + Kies sor ee— 1) 


STEADY STATE ERROR SPECIFICATION 


ERROR SPECIFIED (ERROR J SENPUT) 1S See ee 


BANOWIDTH SPECIFICATION 


NO BANDWIDTH REQUIREMENT IS SPECIFILD. 


meio SPEC I F eR 
CLOSEC LOR SOLU RON 
Aaah oD CCE +C OMEGA = OsIGO sha 08 
NUMBER GF ROCTS = 2 
ROOT NC. ] 2 
Ries —s3F+° ] © See) 


IMAG PT. ef2dt+ 1. =eSi zeae 
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SOLUTION OF THE PROBLEM 


H#ReRKSKEKESEKHEKEKHEEE HEHE FF 


PROGRAM REQUIRES 11 ROOTS FOR SOLUTION. 


STEADY STATE ERROR SOLUTION 


STEADY STATE ERROR IS -1LQ00E+400 


RELATIVE SPECIFICATION 
ERROR IS 3.727E-06 PERCENT 


ERROR SPECIFICATION IS SATISFIED. 


SOLUTION GF THE SYSTEM OF LINEAR LCQUATIONS 


VAR TABLE PARAMETERS RELATIVE ERROR 

VARTABLE ( 1) CD(t) -00 SE OL IP ERCENT 
VARTABLE (  _) CO(1) -1.469%F-09 PERCENT 
VARTABLE ( 3) CD(2) -3.589E-O09 PERCENT 
VARTABLE ( 4) CD( 3) 1.759E=69 PERCENT 
VARTABLE ( 5) CO(4) «2? €-Ol PERCENG 
VARTABLE ( 6) CD(5) sVO>E=0 IPE Gehl 
VARTABLE ( 7) KT#KX#CN (9) =5.242 baOo PERCENT 
VARTABLE ( 8) KT#KX#CN (1) —-l1.2/f “E=95 PERCENT 
VARTABLE ( 9) KT#KX#CN (2) -3.43i1&=08 PERCENT 
VARTABLE (1.4) K T#KX#CN (3) 9.962F=08 PERCENT 
VARTABLE (11) KT#KX#CN (4) F.7T2IE-O7 PERCENT 
VARTABLE (12) KT#KX#CN (5) 3.143£-06 PERCENT 
VARTABLE (13) KT#CN( ©) 2.72 7E-06 BEX CENT 


NOTE - RELATIVE FRROR STATED ABOVE 
MEASURES THE ACCURACY OF THE 
SOLUTION FOR EACH VARTABLE. 
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SOLUTIUN FOR INDIVIOUAL PARAMETERS 


NUMERATOR CF CASCADE COMPENSATION 
CN(S)= 1.20990°9&+0C SC 5) 

+ 6.7668E+01 SC 4) 

+ 1.7/7 (SE Haeees 3} 

+ 2.2639E+04 S( 2) 

+ 1.3997F6+05 S( 1) 


+ 263002 E +05 


ZEROS GF THE FILTER 
ROOT NO. i 2 3 4 5 
REALE PT. —-«2Z0C+ 72 =e l2GtO2 yure lL BE+02 8-. SSE —-.9ZErot 


IMAG PTI. CGEt+20 oe JOE+29 © 00E+900 © J0E+90 e0O0E+U0 


NOTE - ACCURACY OF ROOT FINDING SUBROUTINE 
HAS BEEN CHECKED BY EXPANDING THESE 
ROOTS, AND COMPARING THE COEFFICIENTS 
OBTAINED WITH THE COEFFICIENTS OF THE 
ORIGINAL POLYNOMIAL. 
ERROR NOTED IS LESS THAN ONE PERCENT 
FOR ALL COEFFICIENTS. SOLUTION FOR 
ROOTS IS CONSIDERED ACCURATE. 


NOTE 


CENOMINATOR OF CASCADE COMPENSATIO’: 


CD(S)= L.89G0E+09 St 5) 


+ 7Te2'oDF+O]L SC 4) 


t2e 9 Le+ 3S (es) 


$2. 13512344 Sle) 


t+ LeJi24E4+95 SC 


+ 4216341455 


POLES OF THE FILTER 


ROOT NO. ] 2 3 4 


REAL Pe —e24E+ 32 —~e24E+U2 --89E+tU1] -./8E+91 


IMAG PT. Se, 2 elLLE+H2 eoGe+ -~.26E4+9])] 


- ACCURACY FAS BEEN CHECKED,AS BEFORE. 
ERROR NOTED ITS LESS THAN CNE PERCENT 
FOR ALL CCEFFICIENTS. SOLUTION FOR 
ROOTS ITS CUNSTDERED ACCURATE. 


FORWARD GAIN (KT) = 1.6020E+03 


FEEDBACK GAIN (KX) = 1.2911F+4+99 
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e26E421] 


CLOSED LOOP TRANSFER FUNCTION 


CF SOLVED PROBLEM 


NUMERATOR POLYNOMTAL 


GN(S)= 


ZEROS OF THE TRANSFER 


ROOT 


REAL 


IMAG 


ROOT 


REAL 


IMAG 


Ni). 


PT. 


PT. 


NO. 


PT. 


PT. 


L.6092GE403 


1-6447F4C5 


1.2280t+06 


le 134ee. 6 


22ef652B8F 409 


2e4737t+10 


Le41l2i+11 


4.3966F+]1 


5-8275t+11 


1 


=F 2UEF92 


oCOF+NN 


6 


-.846+0)] 


eO0GE+00 


S( 


S ( 


S( 


S( 


S ( 


S( 


S( 


ai 


3) 


) 


6) 


>) 


4) 


3) 


2) 


1) 


FUNCTION 


2 


—e 1 BE+02 


e46E-32 


7 


e92E+0)] 


eOOE+C 


4 4 
“2 1 GER eee = ol ZEFo2 


~e50F+91 


eG0E+G60 


NOTE - ACCURACY HAS BEEN CHECKED,AS BEFORE. 
ERROR NOTFD YS LESS THAN ONE PERCENT 


FOR ALL COEFFICIENTS. 


SOLUTTON FOR 


ROOTS IS CONSTDERED ACCURATE. 
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5 


-.12E+02 


e 36E-02 


CENOMINATOR 


GO(S)= 


BOE Ses Oir ati tebe Sy Sul EM 
CF THE CHARACTFRISTIC EQUATICN 


ROOT 


REAL 


IMAG 


ROOT 


REAL 


IMAG 


ROOT 


REAL 


IMAG 


NO. 


PT. 


PT. 


NO. 


PT. 


PT. 


NO. 


PT. 


PT. 


PULYNOW TAL 


beeciGes tt eee S{alels) 


lee D5 Ge +02 


Bae Sleigh tees 


2-89C7TE4+%5 


Go13l 8ct56 


lier to 


Le22(Sk+C9 


Bie S4514C9 


5 Alelle ie: sia ae 


2e22G  s+4+1)] 


6.4/50r4+11 


-~,., 29640] 


eo2e+ 1 


11 


are HEC 


eOUEFCS 


S( 


S ( 


S ( 


Sif 


oA 


S ( 


S { 


S{ 


S ( 


AND ROOTS 


—.ouer | 


-~ 5264351 


~~ LSE 02 


10) 


9) 


8 ) 


i. 


6) 


5) 


4) 


3) 


2) 


1) 


052E+09 


2 


7 


SI 


— . 90m, 1 


e52E4C9 


8 


—~e52E4°O 


4& 
-.52E4+90 
9 


-~.18E&+02 


e52E+90 


5 


alee 


e 52E+00 


* 
Ll. 


~ 18E+C2 


~.52E+00 


Example 5,3 
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4 3 2 iw 4 
s + 30s” + 273s” + 740s + 576 O(s) 





Figure 5,3 
Block Diagram for Example 5.3 


The gain of the open loop device shown in Figure 5.3 is less 
than -20 db for all frequencies. Design a filter compensator and 
find values of forward and feedback gain which will provide positive 
gain at low frequencies and a bandwidth of zero db at 30 radians per 
second, 


Dominant roots are specified by € = 0.5 and a LO. 


No steady-state error requirement is specified. 
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LO (ie ae ra 


HSRC HH HEHE 


OPER Hite TOANSFER FUNC T 10 tees) 


NUMERATGR CF iivrN LOOP 


TRANSFER SEU ith = TNS} 


TNGS).=" “62 eet I 


CEND tT ists iP he! LOOP 
TRANSE eeu ie = TOUS} 


Os )= le et 4) 


tesee eet 1 oS (kf) 


ROG er TCS) 


— 
i) 
Ww 
f 


ROUT NU. 


REAL PIs ~el2rtie = 6 LS ae eel oe et we 


IMAG PT. 2° Cc4+22 2256435 SCOEFOG  .C5E+90 
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FEEDBACK PATH TRANSFER FUNCTION - H(S) 


(0) ORDER FEEDBACK IS SPECIFIED. 
NUMERATOR OF FEFDBACK PATH 
TRANSFER FUNCTION - HN(S) 

KN(S)= 1.00C0E+00 

DENOMINATOR OF FEEDBACK PATH 
TRANSFER FUNCTION - HD(S) 

HD(S)= 1.COCOE+00 S( 1) 


+ 5.COGOE+CO 


CLOSED LULIP TRANSFER FUNCTION 
G(S) - UNCOMPENSATED 


NUMERATOR CF CLOSED 0G 


TRANSFER FUNCTION - GN(S) 


GN(S)= KT#*( 6.370 sF+G1 Sf 1) 


+ 3. 1503ErO2 | 


DENOMINATOR CF CLOSED LOOP 
TRANSFER FUNET TON — GOTRS 


GPGS)= 1. eee) 


+ 2.5°. _F+6] S ( 4) 


+ 47 4 ee le 


+ 2. Pear wee 5 (I 


ME ee ey, 


+ ( 2.050643 + KUeRKeG es eee 
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SPEeOY STATE cKROL SPEC PRICAT hs 


NO STEADY STATE ERROR REQUIREMENT 
IS SPECIFIED. 


BANOWIDTH SPECIFICATION 


GAIN OF eCOOE-O1 OB AT FREQUENCY 
OF 3.900E+%!1 RAD/SEC ITS SPECIFIED. 


ROOTS SPECIFIED FOR 
CLOSED LOUP SOLUTION 
ZETA = e5V9CE+00 OMEGA = 1.90909°7E+91 
NUMBER OF ROCTS = 2 
ROOT NOQ. l 2 
REAL PT. -.50E+01 -.50E+491 


IMAG PT. ~87E+0L -.87E+01 


SOLUTION OF THE PROBLEM 


HTHEEEHHEHHEHHEEHREEKEEH REE HE EE 


PROGRAM RECUIRIES 9 ROOTS FOR SOLUTION. 


BANDWIDTH SCLUTION 


GAIN IS eilct-O1 DB AT FREQUENCY 3. ‘O0E+01 RAD/SEC 


Reba IVE EeROR 1S ~JO0E=)!1 “PERCENT 


PANDIV: TOT SPECIFICATION IS SATISFIEO. 
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VAR pA 


--_ 
= 
‘ 


VAK TABLE 


VAR TABLE 


VARIABLE 


VAR TABLE 


VARIABLE 


VAR TABLE 


VARTABLE 


VAR TABLF 


VARIABLE 


tei) 


le SYSTEM? OF 


PARAMETERS 


CO( 0) 


CD(1) 


CD(2) 


CO(3) 


CD(4) 


KT*#KX#CN(%) 


KT#KX#CN (1) 


KT#KX#CN (2) 


KT#*KX#CN (3) 


KT#KX#CN (4) 


NOTE - RELATIVE ERROR STATED ABOVE 
MEASURES SHE ACCURACY OF THE 


SOLUT 10% 


FOR EACH VARIABLE. 
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LINtr AR EQUATIONS 


RFLATIVE ERROR 


eo James 


— jee O2E=Oo 


= fie DiGtl E=O8 


= Sie LE Oy, 


1.559E-06 


7.654E-05 


leceeesk 04 


Ge 762. —035 


5223 IE-02 


2 -441E-Ol 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


SCLUTTCN FOR THDIVIDUAL PARAMETERS 


NUMERATOR uF CASCADE COMPENSATION 


CNUS T= | ong. S04) 


+. Bey 7 ela ee 


+ 5.47141 $22 Stee 


+ 3.497345) ee SC) 


+ 0). 4 Oe 


PERS leit 0 el fea 


ROOT Nt. ] 2 4 4 


REAL PT, —~sliZzer <2 - 1SE ez - 61649) =. OlE+ a 


IMAG Peli-< ai ee < eo) JE 1% -.-36L+)1 ao 5 tea 


NOTE = ACCURACY (F "(4 ] FIADIN Sur ieee le 
H4S BEEN CHECK eb BY TeX Pe De a oe 
REOTS, VID CGMP AR TN eliGlge ce Oteeteelee = Nias 
OBTALNED wlth “(ht Cai Beer Se eee 
ORT eli Ak Put Y NOM 1 Alle 
EMRIR NOTED 1S LESS# VRANE CRE PERCENT 
FOR Alb COCLEF LCT CN Tie Sie eee ie Gis mete 
ROWTS [Ts CUcSIDERED ACEUA a TE. 


59 


CENOMINATOR OF CASCADE COMPENSATION 


CO(S)I= 1.0924t+00 St 4) 


+ 1.25C0E+02 S( 3) 


toe 4ot+93 Sl. 2) 


4 1.6789E+95 S( 1) 


+ 264427E+06 


pees OF THE FILTER 
ROOT NO. 1 2 B 4 
BeMiePpt. --476¢+02 =.4782 -.15t4+02 -.1er Zz 
[MAG@PT = = litt 2 - Whee —-. 27+ t 02 s2reT og 
NOW eae OURACY HAS BEEN CHECKED,AS BEFORE. 

ERROR NCTED TS LESS THAN ONE PERCENT 

FOR ALL COEFFICIENTS. SOLUTION FOR 

ROOTS IS CONSIDERED ACCURATE. 
FORWARD GAIN (KT) = 4-2411E+04 


FEEDBACK GAIN (KX) = 7.1267E+906 
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NOTE 


CLUSED LOUP TRANSFER EUG iia 


(F SOUVED PREBG@Ee 


NUMERATOR POLYROMEAL 


CNGS3= 2267 Vee ast 


+ Lebhigot 4 raat 


+ 1. CI 9S 


+ 1.53 46r4 S ( 


+ 6% 55.3] | Gal eS 


Pda ae) 


5) 


4 ) 


3) 


2) 


1) 


ZEROS (.F TH= TRANSFER FUNGCTETE 


RGOT Wu. ] 


RE APT. 


[MAG PI. -.36c+4° } 


Tee ACY HASMBEEN 


= 5 Olea —. SlE7 


De} 


. 3 GE tel oO 


CHECKED, AS HERDRig. 
ERROR NCTED FS GREATER THAN ONE PERCENT 


Gk dle LFASTOONL COEPRICTERT. 


POR RCOTS IS SOT CONSIDEREG BCCUR AGE. 
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oad Sie 4 7 


4 > 


elie Bites Diet 


GUEtOC eUDEFS0 


CENOMINATOR POLYNOMIAL 
GO(S)= 1.C7CGE+00 S( 9) 
+ 1.60C0E+02 S( 8) 
teed i oZte ros SCT) 
+ 4.4527E+U5 SC 6) 
+ 1.1274E+07 SC 5) 
+ 1.8947E+98 S( 4) 
+ 2.1474E+09 S( 3) 
+ 1.62C1F4#190 Sf 2) 
+ 7.5736F+10 Sit ole) 
+ 1.6983E+11 
POLES OF THE SYSTEM AND ROOTS 
CF THE CRARACTERISTIC EQUATION 
ROOT NG. 1 2 3 4 5 
REAL PT. -.5GE+01 -.50E€+0T -.15E+02 -.156+02 -.20E+02 


IMAG PT. sO7TE+O] “=STETEFTL ~87E+00 -.87E+90 e87E+C0 


ROOT NO. 6 qT 8 9 
REAL PT» --e20E+O02 -.25E&+02 -.256€+02 -.30E+092 


IMAG PT. —--@7E+0C o87E+90 -.-87E+900 -00€+00 
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INPUT DATA 


SCHERER HEHEHE EE 


OPEN LOCP TRANSFER FUNCTION -— T(S) 
NUMERATOR OF OPEN LOOP 
TRANSFERS FUNCTIONS - @N(S) 
NUMBER OF ROCTS = 2 
ROOT NC. l | 2 
REAL PI. --.60L+00 €.13E+ 2 


IMAG PT. eal + Oe et 13 


COERFFICTERTS @P GAP ANDE Bm PUL YNOMIAL 
TN(S)= lec eetrao Si 24 
tela Viet 1 SE de 


+ F890 9F49] 
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CENOMINATOR (FF OPEN LOOP 
TRANSFER EUR 1 Cli eDitess) 


NUMBER. OR SRUG 1S =e 


ROOT Ni. l 2 3 4 


REAL PT. lobe =. SIR) sae oe te See tees 


IMAG PT. o! eet SCE). =. 30! Wal oo OEE ow 


ROOT NC. A 


ie et cles «1 i td 


TMAG PT 6 rh Re) 


COErrFiC Vent S Uf’ x2aNheo Pol YNOMiAG 


BIOS = eG Eee 5S (6) 


ee Oe) RG oy 


eo Ne 


+> Geol y ote oO ote} 


fede rour. 4 st 2) 


+ 8.°%°4345474 S( 1) 


frre el oe C=) 


65 


—e1l5rt+u2 


~80E+C] 


FEEDBACK PATH TRANSFER FUNCTION - H(S) 
(0) ORDER FEEDBACK IS SPECIFIED. 
NUMERATOR OF FEEDBACK PATH 
TRANSFER FUNCTION - HN(S) 
HN(S)= 1.-9090C0E+00 
DENOMINATOR OF FEEDBACK PATH 
TRANSFER FUNCTION - HO{(S) 


HD(S)= 1.9000E+00 
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CLOSEC LOGP TRANSFER FUNCTION 
G(S) - UNCCMPENSATED 


NUMERATOR (CF CLOSED LOOP 

TRANSFER FUNCTION - GN(S) 

GN(S)= KT«#{ lL. whJIEede St 2} 
+ | <GU@me+oy S({ 1) 


PENOMATOR CF CLOSED LBOr 
TeemoteR FUNCTICNS=] GO{ S$} 


CDiSd=. lect Jeeves 6.) 
# 4.9 . ott] St 5) 
+ 9.35051 +02 SC 4) 
+ GS.5170) 403 SE 3) 
+ ( 4.5 168424 + KT#KX# 1.00 SEF o) SC 2) 
+ ( &8.8424664.4 + KT#kKX#® 1.965084 °1) S( 1) 


AE ( wc — = KT#KX # 7.~8QGQLGEF '1) 
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STEADY STATE ERROR SPECIFICATION 


ERROR SPECIFIED (ERROR / INPUT) IS .LCGCE+00 


BANCWIDTH SPECIFICATION 


GAIN OF -3.0(C0E+0G0 DB AT FREQUENCY 
OF 2.COCE+(1 RAD/SEC IS SPECIFIED. 


ROOTS SPECIFIED FOR 
CLOSED LOUP SOLUTION 
ZETA = .797CE+C0 OMEGA = 4&.0000E+00 
NUMBER OF ROOTS = 2 
ROOT NO. 1 Za 
REAL PT. -.-28E+01 -.28E+01 


IMAG PT. 928EtG)] . ~s2 SE L 
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SQLUTICN UF THE PROBLEM 


HHEKHREHEHEKEKHEHEEEEE HEHE EH 


PROGRAM REQUIRES 11 ROOTS FOR SOLUTION. 


BANOWIDTH SOLUTION 


ITERATIVE PRCCEDURE FOR BANDWIDTH SOLUTION 


PROCESS IS CONVERGING 


ITERATION GATN-DB RELATIVE &RROR 
31 — 26 (ae to) 1. 295 MEISE ERC BINg 
36 <6 ee Oo 6% 6G SE— J 1 Pine ent 


GALN SIS —3.0206+''0 DOB ATIeEREQUENCY 2. GUEt0 LUNAR AGEs 


Bel ATVEwERROR IS -6.6556-9] PERCENT 


BANDWIDTR SPECIFICATION IS SATISFIED. 


STEADY STATE ERROR SOLUTION 


STEADY STATE ERROR IS _.1OGm Bate 


RELATIVE SPECIFICATION 
ERR SIS 3.274F-08 PERCENT 


PreeReoPECTEICATION 1S SATISEFED. 
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SOLUTION OF THE SYSTEM OF 


VAR TABLE 


VARIABLE 


VAR TABLE 


VAR TABLE 


VARIABLE 


VARTABLE 


VAR TABLE 


VARTABLE 


VARIABLE 


VAR TABLE 


VAR TABLE 


VARTABLE 


VARIABLE 


VARTABLE 


( 


1) 


2) 


a 


4) 


5) 


6) 


1) 


8) 


9) 


(19) 


(11) 


(12) 


(13) 


PARAMETERS 


CD(0) 


CO(1) 


CD(2) 


CO(3) 


CD(4) 


CO(5) 


KT#KX#CN (0) 


KT#KX#CN(1) 


KT#KX#CN (2) 


KT#KX#CN (3) 


KT#KX#CN (4) 


KT#KX#CN(5) 


KT#CN( 9) 


NOTE - RELATIVE ERROR STATED ABOVE 
MEASURES THE ACCURACY OF THE 
SOLUTION FOR EACH VARIABLE. 
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LINEAR FQUATIONS 


RELATIVE ERROR 


=2.2979R- 069 
eQQIE-01 
eJJDE-O1 
1.689E-09 

eGOIJE-O1l 
3.593E=00 
1.932E-08 
EPG HE) my 0) 
1.55 5E=J6 
1.347E-05 
8.9095E-05 
4 2292E-04% 


=—5. 1 Zee +0% 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


PERCENT 


NOTE 


SCLUTION FOR INDIVIDUAL PARAMETERS 


NUMERATOR CF CASCADE COMPENSATION 


CN(S)= so OPPs tao 


“+ Ge.57TSGE FS 4 


+ 8.3424 +52 S( 3) 


4 6.9979, +53 St 2) 


492 (ONDE + 4a. SA a) 


Zee sme Ur 14%. FULT ake 


ROOT NO. ) 4 


N 
lu 


REAL Pl. -elorame -. SETS =—290 et ae oc 


Meee fee —. fOr +s | o Cerra eNOL4CO ——s27ere 


= ACCURACY Che RUD teRINDENG SMRROUIeLME 
HAS BEF: CHLCKED BY EXPAN DIGeyTHeSE 
ROOTS, AND COSPARING THE CCEFFICIENTS 
OBTAINED FITH THE CORRE ite sur ine 
ORIGINAL FPOLYNOMTIAL. 
ERROR NOTED IS LESS THAN CNE PERCENT 
POR ALL COEFFICIENTS. SHEL UTION + Oe 
ROOTS IS CONSIDERED ACCURATE. 
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~.-34E+0] 


sere rut 


CENOMINATOR CF CASCADE COMPENSATION 
CO(CS)= LeNVQO0E+NO SC 5) 

+ 1.7057E+02 St 4) 

+ 1.1647E+04 S( 3) 

+ 4, COOTE+O5 oS (eZ) 

+ 4.)9045E+96 S( 1) 


+ 5.85 76E+06 


POLES OF THF FILTER 
ROOT NO. 1 2 3 4 5 
REAL PT. -.806+02 -.38E+02 -.38E+02 ~-.18Et+01 -.13£+02 


IMAG PT. eGCGOE+O00 -—-.43E+02 © 43E+02 eOGE+00 ~O0E+00 


NOTE ~ ACCURACY HAS BEEN CHECKED,AS BEFORE. 
ERROR NCTED IS LESS THAN ONE PERCENT 
FOR ALL COEFFICIENTS. SOLUTION FOR 
ROOTS IS CONSIDERED ACCURATE. 


FORWARD GAIN (KT) = 2-6980E+C6 


FEEDBACK GAIN (KX) = 1.21111E+00 


a 


NOTE 


CLOSED LOCP TRANSFER FUNCTION 
CF SOLVED PROBLEM 

NUMERATOR POLYNOMIAL 
GN(S)= 2.6980i+74 St 7) 

+ 1.7464'498 S(t 6) 

Ae oe Oe 9D Oo 

A Pere Le Sl 6G) 

+ Ge gosr* ll! St. 3) 

Seen cele 1 2S 2} 

tee 24 Oe ee S (1) 


# 1.0°5914+13 


Coe SINE TRAN OF cR + URerren 
RCOT NQ. l 2 3 4 


REAL PT. =e lout oO -~.LoEt32 ~el5E+C2 -~ 69E+ Ol 


[IMAG PT. .C5L405 -.76E451 =. TOE4O1 «=. COE FDC 


ROOT NG. 6 q 

REAL ol ees —~ - 324c+G]) ~  34E+uU1 

PMAGePT. ~edfit.it e276 +)1 

- ACCURACY HAS KEEN CHECKED, AS BEFORE. 
ERRUR NCTED IS LESS THAN CNE PERCENT 


FOReARY CCEFFIGRENTS. SOWMUTTON FOR 
ROUTS IS CONSIDERED ACCURATE. 
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~90E+00 


CENOMINATOR POLYNOMIAL 


GDO(S)= 


POLES OF THRE SYSTEM AND ROOTS 
CF THE CRARACTERISTIC EQUATION 


ROOT 


REAL 


IMAG 


ROOT 


REAL 


IMAG 


ROOT 


REAL 


IMAG 


NO. 


PT. 


PT. 


NO. 


PT. 


PT. 


NO. 


PT. 


PT. 


1.960C0E+00 S(11) 


2-1957E+02 


2-CIBBEt+O4 


1.1489E+06 


3e-9/23E+07 


9.-08COE+O08 


1.3919F+10 


1-4210E+11 


9-4/53t411 


3-9646E+12 


Ge-6743E+12 


ieee) 7 tls 


1 


~28E+91 


6 


ll 


as ~43E+02 


e00E +90 


$(10) 


S { 


S{ 


S ( 


S ( 


Sf 


S({ 


S 


S ( 


S{ 


- .28E+01 


9) 


8 ) 


7) 


6) 


2) 


4) 


3) 


2) 


1) 


e28E+09 


2 


7 
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3 


=e 85E+01 


e28E+09 


8 


~.26E+92 


—~e2oe wm 


4 
-.85E+01 


a) 28E45% 


9 
~.34E+92 


e28E+09 


5 


-. life 


e28E+tO0 


~, 


i 


—.34E+02 


6, Conclusion 

The method presented has been shown to be fast and accurate, 
and capable of designing systems which will meet specifications 
for root dominance, steady-state error, and bandwidth. 

The principal disadvantage of the method presented, exclusive . 
of the restricticos mentioned in Section 4, is the large number of 
filter sections required, particularly for higher order open loop 
systems, The filter sections tend to further complicate the 
originally complex systems, 

The possibility of reducing the number of filter sections 
employed in the final design is worthy of further study. Approxi- 
mation of the characteristics of the cascade compensation may be 
accomplished with a filter of lower order, but analysis of the 
effects of such approximation on system specifications is beyond 


the scope of this paper. 
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